function obj = coef_J_pure_pTop_func_new(in1,in2,in3)
b1 = in2(1,:);
b2 = in2(2,:);
b3 = in2(3,:);
n1 = in3(1,:);
n2 = in3(2,:);
n3 = in3(3,:);
r1 = in1(1,:);
r2 = in1(2,:);
r3 = in1(3,:);
t2 = b1.*n1;
t3 = b2.*n2;
t4 = b3.*n3;
t5 = n1.*r1;
t6 = n2.*r2;
t7 = n3.*r3;
t8 = n1.*r2.*2.0;
t9 = n2.*r1.*2.0;
t10 = n1.*r3.*2.0;
t11 = n3.*r1.*2.0;
t12 = n2.*r3.*2.0;
t13 = n3.*r2.*2.0;
t14 = -t5;
t15 = -t9;
t16 = -t6;
t17 = -t11;
t18 = -t13;
t19 = -t7;
t20 = t8+t9;
t21 = t10+t11;
t22 = t12+t13;
t23 = t5+t6+t7;
t27 = t2+t3+t4;
t24 = t8+t15;
t25 = t10+t17;
t26 = t12+t18;
t28 = t5+t6+t19;
t29 = t5+t7+t16;
t30 = t6+t7+t14;
obj = [t23.^2,t23.*t26.*2.0,t23.*t25.*-2.0,t23.*t24.*2.0,t23.*t30.*-2.0+t26.^2,t20.*t23.*2.0-t25.*t26.*2.0,t21.*t23.*2.0+t24.*t26.*2.0,t23.*t29.*-2.0+t25.^2,t22.*t23.*2.0-t24.*t25.*2.0,t23.*t28.*-2.0+t24.^2,n1.*t23.*2.0,n2.*t23.*2.0,n3.*t23.*2.0,t23.*t27.*-2.0,t26.*t30.*-2.0,t20.*t26.*2.0+t25.*t30.*2.0,t21.*t26.*2.0-t24.*t30.*2.0,t20.*t25.*-2.0-t26.*t29.*2.0,t20.*t24.*2.0-t21.*t25.*2.0+t22.*t26.*2.0,t21.*t24.*2.0-t26.*t28.*2.0,n1.*t26.*2.0,n2.*t26.*2.0,n3.*t26.*2.0,t26.*t27.*-2.0,t25.*t29.*2.0,t22.*t25.*-2.0-t24.*t29.*2.0,t22.*t24.*2.0+t25.*t28.*2.0,n1.*t25.*-2.0,n2.*t25.*-2.0,n3.*t25.*-2.0,t25.*t27.*2.0,t24.*t28.*-2.0,n1.*t24.*2.0,n2.*t24.*2.0,n3.*t24.*2.0,t24.*t27.*-2.0,t30.^2,t20.*t30.*-2.0,t21.*t30.*-2.0,t29.*t30.*2.0+t20.^2,t20.*t21.*2.0-t22.*t30.*2.0,t28.*t30.*2.0+t21.^2,n1.*t30.*-2.0,n2.*t30.*-2.0,n3.*t30.*-2.0,t27.*t30.*2.0,t20.*t29.*-2.0,t20.*t22.*2.0-t21.*t29.*2.0,t21.*t22.*2.0-t20.*t28.*2.0,n1.*t20.*2.0,n2.*t20.*2.0,n3.*t20.*2.0,t20.*t27.*-2.0,t21.*t28.*-2.0,n1.*t21.*2.0,n2.*t21.*2.0,n3.*t21.*2.0,t21.*t27.*-2.0,t29.^2,t22.*t29.*-2.0,t28.*t29.*2.0+t22.^2,n1.*t29.*-2.0,n2.*t29.*-2.0,n3.*t29.*-2.0,t27.*t29.*2.0,t22.*t28.*-2.0,n1.*t22.*2.0,n2.*t22.*2.0,n3.*t22.*2.0,t22.*t27.*-2.0,t28.^2,n1.*t28.*-2.0,n2.*t28.*-2.0,n3.*t28.*-2.0,t27.*t28.*2.0,n1.^2,n1.*n2.*2.0,n1.*n3.*2.0,n1.*t27.*-2.0,n2.^2,n2.*n3.*2.0,n2.*t27.*-2.0,n3.^2,n3.*t27.*-2.0,t27.^2];
